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Abstract 

In this paper based on the notion of entanglement witness, a new measure of entan- 
glement called floating entanglement witness measure is introduced which satisfies some 
of the usual properties of a good entanglement measure. By exploiting genetic algo- 
rithm, we introduce a classical algorithm that computes floating entanglement witness 
measure. This algorithm also provides a method for finding entanglement witness for a 
given entangled state. 

AMS: 65C99; 81P99 
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Entanglement is one of the most interesting properties of quantum mechanics and the key 
resource of some quantum information and quantum computation processes, such as teleporta- 
tion, dense coding and quantum key distribution [T, ^2] . However characterizing and measuring 
the entanglement has tantalized physicists since the earliest days of quantum mechanics, and 
even today there is no general qualitative and quantitative theory of entanglement. 

Among the known criterions for distinguishing between separable and entangled states 
[3 m m El 13 [HI E], entanglement witnesses (EWs) have a special importance for detecting 
the presence of entanglement. Finding EW that leads to solving separability problem, is an 
interesting but computationally a demanding job as it has been shown that the separability 
problem lies in the class of NP-hard problems [10] . There has been a lot of efforts for finding 
EW for given quantum state; some of them propose EW for some special cases [TT] while 
others try to find approximate solutions by using different methods like semidefinite and linear 
programming [121 [IS] • In addition to the question of " Is it entangled?" there is the question 
of "How much entangled?". Some measures of entanglement related to EW are introduced by 
Bertlemann et al [H] , Brandao and Vianna [T5l [16] . 

The aim of this paper is two fold: firstly by defining a slightly different definition than 
of EW which is called floating entanglement witness (FEW), we introduce a new measure of 
entanglement namely FEW measure. Secondly based on the genetic algorithm (GA) we offer 
an algorithm which computes FEW measure and also finds EW for every entangle state. GA is 
a powerful and intelligent technique for global optimization, adaptation and search problems 
[T7] [T8| [To] . It reveals its power especially when we are dealing with highly nonlinear and 
large search spaces. This technique is inspired by natural evolution of species which is based 
on selection, inheritance and mutation. Readers are referred to appendix I for an overview of 
GA. 
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The paper is organized as follows: In section 2 we briefly recall the definition of the EW 
and define our concept of the FEW measure. In section 3 an algorithm is proposed for finding 
FEW measure and EW for a given entangled state. Section 4 is devoted to some examples 
such as: Bell, Werner, mixture of GHZ and W and one parameter two-qutrit states. The paper 
is ended with a brief conclusion and two appendices. 



As mentioned in the introduction one of the methods for detecting entanglement is applying 
EWs. Let us first recall the definition of entanglement, separability and EWs. A density matrix 
p is called separable or unentangled if there are positive pi's with 'YliiPi = 1 ^'^^ product states 



otherwise it is called entangled. 

Definition 1. A Hermitian operator W G B{l-Ld^®'Hd2 ■ ■ ■®'Hd^) (the Hilbert space of bounded 
operators) is called an EW detecting the entangled state Pe if Tr{yVpe) < and Tr{yVps) > 
for all separable state ps G S. 

Therefore, if for state p we measure TriyVp) < 0, we can be sure that p is entangled. 
This definition has a clear geometrical meaning. Thus, the set of states for which Tr(VVp) = 0, 
is a hyperplane in the set of all states, that cuts this set into two parts. In the part in which 
TriyVp) > 0, hes the set of all separable states. The other part (with Tr(>Vp) < 0) is the set 
of entangled states detectable by W. From this geometrical interpretation it follows that for 
each entangled state Pe, there exist an entanglement witness detecting it. This statement is 
proved in [4]. 



2 FEW measure 



(^i ) such that 




(2.1) 
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Definition 2. A traceless Hermitian operator Z is called a floating entanglement witness 
(FEW) for detecting the entangled state pe if Tr{Zpe) < minp^g^ Tr{Zps) where S is the set 
of separable states. 

It is necessary to note that the existence of Z comes from the existence of EW for every 
Pe- Due to the convexity of the separable states, the minimum in the above definition comes 
from the border of separable states (pure product states). One can easily show that 



is an EW since it satisfies botli conditions of EW, i.e.: Tr(VVps) = Tr{Zps) — p > for all 
separable ps and Tr(yVpJ = Tr{ZpJ —p < Tr{Zps) —p < 0, which means that Tr(>VpJ < 0. 

A closely related problem to the EWs is characterization or quantification of entanglement 
by EWs (see [16]). Based on the concept of FEW we introduce a new computable FEW 
measure for quantifying entanglement of a given quantum state: 



W = Z-pI 



where p = mmTr{Zps), 



(2.2) 



E{p) := max < , max 




where 



0, 11^112 = 1}, 



and II II 2 := \/Tr{Z'^ Z\ denotes the Hilbert-Schmidt norm, 
following usual requirements of an entanglement measure pO] 



The FEW measure fulfills the 



Proposition 3. -E'(p) satisfies the following properties: 



(i) For every separable state o"s G 



S, E{as) = 0. 



(ii) Local unitary operations leave E{p) invariant, i.e. 



E{p) = E{Ul (g)...(g) UlpUi ® . 



(2.4) 
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(iii) FEW measure dose not increase under local operations and classical communication 
(LOCC) protocols [21], i.e., 

E{£ip)) < E{p). (2.5) 

(iv) The FEW measure is a convex function, i.e., 

E{\p + (1 - \)a)) < XE{p) + (1 - X)E{a). (2.6) 

(v) FEW measure is a continuous function. 

Proof: (i) It is easy to see that for every Z we have miup^ Tr{Zps) — Tr{Zas) < hence 
Eias) = 0. 

(ii) (12. 4p follows from the invariance of S and ||.||2 under local unitary operations. 

(iii) Although we did not find rigorous proof for monotonicity of FEW measure under general 
LOCC maps S, it can be proved for isometry ones for which we require S'^S = ££'^ = I. To 
see this we note that 

E{8{p)) = max |o , max [mmTr{Zps) — Tr(^Z£{p)^'\ j- (by definition) 

= max < , max [ min Tr(^£\Z)£\ps))—Tr(^£\Z)p^~^ > (isometry) 

= max jo , max^^ [mmTr{Z'p,) - Tr{Z'p)] | {£K^) 

< max 1 , max \mmTr{Zps) - Tr{Zp)] I = E{p) i£\A) C A) 
zeA pses ■' J 

where in the second equality we use the fact that £"^(3) = S, [22] and Tr (^Z£(p)) = Tr (^£'^{Z)p) . 

(iv) To prove the convexity we note that 

E{Xp + (1 - X)cr)) = max i , max [minrr(Zp,) - Tr{Z{Xp + (1 - X)a)] \ 

I Z^A Ps^S \ 

= max < , max \XmmTr(Zps) + (1 — X) iamTr{Zps) — XTr(Zp) — (1 — X)Tr(Za)] 

ZeA ps£S psGS 
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< Amax|o , max [mmTr{Zps) — Tr(Zp)] j> + (l— A) max |o , max [mmTr{Zps) —Tr{Za)] 

where the inequahty comes from the fact that, maxn(/i + f2) < maxfj /i + maxa /i, for every 
function fi, /2 defined on a given set Q. 

(v) To prove continuity, for entangled states p and a we suppose that E{p) > E{a), without 
loss of generality. By definition 

E{p) = max [minrr(Zp,) - Tr{Zp)] = Tr{ZpPs^) - Tr{Zpp), 

where Zp and ps^ are the Z & A and ps & S for which the maximum occurs in the above 
equation. On the other hand we have 

E(a) = max [mmTr(Zas) - Tr(Za)] > VZ [mmTr(Zas) - Tr(Za)] 

=^ E{a) > Tr{ZppJ - Tr{Zpa) 

Therefore 

E{p) - E{a) < -Tr{Zpp) + Tr(Z^a) = Tr(Z^(a - p)) < \\ZpMa - pW^ = \\a - ph = e 

for some real number e > 0. In the last inequality we have used Cauchy-Schwartz inequality 
for Hilbert-Schmidt distance. 

Another usual property of every measure of entanglement is the additivity E[p®^) — nE{p) 
or subadditivity E{p (g) cr) < E{p) + E{a) problem which for FEW measure remains open for 
debate. 

3 Finding FEW measure and EWs using GA 

This section deals with application of GA in finding FEW measure and EWs. At the first 
stage the problem is finding a FEW for a given p. The fitness function is defined as 

JP(p, Z) := min Tr(Zp,) - Tr{Zp). (3.7) 

Pa 
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The task is to maximize JF over Z by using the GA, which imphes that GA tends to make 
Tr{Zp) smaller than mm.Tr{Zps), as well as it increases the difference between these two terms 
in each chromosome. As it is clear from the fitness function we need a second optimization 
procedure that finds min Tr{Zps), for each chromosome. In this subprogram, quasi-Newton 
(QN) optimization method is used by employing optimization function UMINF of IMSL math 
library. For this reason we calculate Tt{psZ) for a large number, iV^jvn '^^ random p^'s then 
use A^Qjv2 of them which have smaller Tt{psZ) as initial points for running QN. The smallest 
resulting value is chosen as niiup^ Tr{Zps). 

Regarding the above considerations the algorithm goes as follows: 

1. Read input density matrix p. 

2. Populate an initial pool of random Z^s (chromosomes), i.e. produce Z^s with random 
parameters. 

3. (a) Find /i := min^^ Tr{psZ) for each Z by using QN. 

(b) Find /2 := Tr{pZ) for every Z in the pool. 

(c) Compute the fitness function T = j\ — /2, for every Z. 

4. Produce the new generation by doing selection, crossover and mutation. 

5. Go to step (IH]) until the stop criteria is met, i.e., maximum number of iterations is run. 

6. (a) Select the chromosome {Z\ with maximum T 

(b) If J" > then T is FEW measure and Z is FEW. 

(c) Compute W = 2 — /iJ (with p = fi) as an EW for detecting p. 

7. if JF < 0, p is separable. 

It is useful to mention that even if we just want to find EW, it is better to find it through 
FEW. An EW, W, must satisfy the conditions in definition [H i.e. a W which Tr(>Vp) < 
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for the given p, subjected to the constraint Tr(Wps) > for all separable states p^'s. The 
latter condition puts our problem in the field of constrained problems. GA is adaptable to the 
constrained problems by defining a proper fitness function. But the constraint of the definition 
[T] slows down the algorithm and therefore instead of it we use definition [2] which breaks the 
constraint and greatly speeds up the algorithm. 

The above algorithm can find a FEW, an EW and FEW measure for any arbitrary entangled 
density matrix. The robustness of the algorithm will be clear in the next section with various 
examples. 

4 Some examples 

In this section we discuss some interesting cases which can help to clarify the subject. In 
all cases we use the fact that any traceless Hermitian operator acting on the Hilbert space 
'Hdi ®'Hd2---®'Hd,^ can be expressed by identity operator Idid2...d„ and generators of Lie algebra 

su{di), Xf'\ as 



where Xq = Id, and Ti^i2...i^ G 

Typically in GAs, there is no strict rule for choosing parameters of them. Probability of 
crossover is usually chosen 0.7 and probability of mutation P^, is a very small number in 
10"'^ order. Number of population N^^ is usually several hundreds to several thousands. 

In the following examples, crossover is two-point, Pc is 0.7, Pm is 0.007 and N^^ is chosen 
20 (for 2 qubit states) and 10 (for other states) times bigger than the number of parameters 
of the problem. Convergence of GA in our problem starts, depending on the quantum state, 
from 20 iteration to higher ones. Therefore the maximum number of generations G^^, in most 
examples is chosen 300 which is big enough for finding proper solutions. Also, there are no 



4-14-1 di-i 




(4.8) 
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lower and upper limits for parameters of QN except that choosing very small numbers make 
QN less successful in finding global minimum and very big numbers make the algorithm very 
slow. In this program, Nqni is chosen about 100 times bigger than the number of parameters of 
the general form of pure separable states and NQi\f2 is chosen about this number of parameters. 

4.1 Two-qubit systems 

It is important both theoretically and experimentally to study entanglement of qubit systems 
and to provide EWs to verify that in a given state, entanglement is really present. 

Here we find an EW for a given two-qubit density matrix. The most general form of a 
traceless Hermitian operator in the space of two-qubit states can be written as: 

3 

Z = ^ Tijai (g) CTj, Too = 0, (4.9) 

i,j=0 

where Tij G M, o"o = I2 and cr^'s with i = 1, 2, 3, are usual Pauli matrices. The equation (14. 9 p 
has 15 parameters (Ty's) which have to be determined by GA. The range of r^j's is determined 
by the condition \\Z\\2 = 1 in every example. Each of the Tj^'s is encoded in a 15 bit binary 
number. Therefore every Z can be encoded in a chromosome of 225 bits. Parameters of 
each chromosome are passed to the QN. In the QN, for finding /i, it is sufficient to take this 
minimum over the p^'s in the form pg = (8> \4'2){4'2\ which can be parameterized as: 

|0,) = cos(/3,)|O) + e*"^sin(/3,)|l), j = l,2 (4.10) 

where aj G [0,27r] and jSj G [0,^]. QN algorithm minimizes Tr{Zps) by determining these 
four parameters. 

4.1.1 Bell states 

We begin with the simplest case that is finding an EW for pure Bell state iV^oo) (V'ool, where 
iV'oo) = 7f (|00) + 1 11))- Crossover is two-point, probability of crossover (Pc) is 0.7, probability 
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of mutation (Pm) is 0.007 and numbers of population and generations for GA are A*", 



GA 



350 



and G^^ = 80 respectively and parameters of QN are Nqj^i 
program gives us the following solution: 
/ 



400 and N, 



QN2 



5. The 



W 



V 





0.296 
0.280 

-0.289 - O.OOli 



0.296 0.280 -0.289 + 0.001i 

0.575 0.288 - O.OOH 0.283 

0.288 + O.OOli 0.578 0.292 



0.283 



0.292 







with E{p) = 0.577. If we don't impose the condition \\Z\\2 = 1, we have: 



\ 





-omsi 

-0.002Z 

-1 - 0.014i 



0.003f 0.002Z 
0.999 

0.999 
-0.002i 



-l + 0.014i 
0.002i 





This is very similar to the EW corresponding to reduction map 

( -1 \ 



Ted 



I - 2|V^oo)(V'oo| 



10 
10 
y -1 / 



Of course choosing any Bell state 



:= o\®oI\^qq), i,j = 0, 1, 



the GA yields an EW similar to reduction EW, Wred = I — 2|V'ji)(V': 



(4.11) 



(4.12) 



4.1.2 Werner states 

One of the most important degraded Bell states is Werner state |24j. A Werner state in 2 ® 2 
system takes the following form: 



1 - F 

Pw = F|^/'oo)(^oo| + (|^io)(^io| + |^oi)(^oi| + 



< F < 1 



(4.13) 
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where \ipij)^s are Bell states defined in Eq. (14.121) . The Werner state is characterized by 
a single real parameter F called fidelity. This quantity measures the overlap of Werner state 
with a Bell state {ipoo)- The is separable for < F < ^ and is entangled for | < F < 1. In 
these cases which the EW is complicated to write down, we give only the FEW measure (see 
Fig. 1). 

4.2 Three-qubit systems 

The procedure for three-qubit states is similar to the one described in the previous subsection. 
It is clear that every three-qubit traceless Hermitian operator Z can be written by tensor 
product of Pauli operators as 

4 

Z = ^ Tijkai (g) aj (g) (jfc, Tooo = (4.14) 

i,j,k=0 

where the number of parameters Tij^ G M, is 63, so every Z can be encoded in a string of 
15 X 63 = 945 bits. Similar to two-qubit case, we work with pure separable states ps = 
|0i)(0i| ® |02)(02| ® |03)(03|, for finding miup^ Tr{Zps), where is defined the same as Eq. 
(I4.10p . Therefore we have 6 parameters for QN. 

For example consider the mixture of |W) = (|100) + |010) + |001))/V3 and |GHZ) = 
(|000) + |111))/V2: 

Pg = g|GHZ)(GHZ| + (1 -g)|W)(W|, 0<g<l (4.15) 

Results for FEW measure which are shown in Fig. 2 are in agreement with witnessed entan- 
glement [16] . 

4.3 One parameter two-qutrit state 

The description of qutrit systems is very similar to the one for qubits. For two-qutrit states 
the operator Z can be expressed a.s Z = Yl^j=o '^ij^i ® Aj, Too = 0, where Aj's with i = 1, ...,8 
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are the well-known Gell-Mann matrices (see appendix II). The number of parameters r^j G M, 
is 80, so every Z could be encoded in a string of 15 x 80 = 1200 bits. We work with pure 
separable states = |0i)(0i| ® |02)(02|, for finding minTr{Zps), where 

= e^''^ siniOj) sin(0j)|O) + e'^^ sin(^j) cos(0j)|l) + cos{9j)\2), j = 1,2 

and r]j,$,j G [0,27r] and 9j,(f)j G [0, |]. As last example consider the one parameter two-qutrit 
density matrix 

P. = + + 2<«<5 (4.16) 

where 

10+) = 73(100) + 111) + |22)) 

a+ = |(|01)(01| + |12)(12| + |20)(20|) (4.17) 

a_ = |(|10)(10| + |21)(21| + |02)(02|) 

It has been shown that is separable if 2 < a < 3 and entangled if 3 < a < 5, FEW 
measure for these states is shown in Fig. 3. 

5 Conclusion 

We have introduced a new measure of entanglement called FEW measure, to quantify entan- 
glement of quantum states which is based on a slightly different definition than of EW. Some 
properties of every measure of entanglement were proved for this measure. Using GA we pro- 
posed an algorithm to compute FEW measure for multipartite systems. For several examples 
including two-qubit, three-qubit and two-qutrit systems, the results of this method were illus- 
trated. As in the examples, this method is capable of finding numerically these operators and 
quantity for every quantum state with any number of dimensions and parties. This algorithm 
also provided a method for finding EW for any entangled state. 
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Appendix I 

GA: GA is a search algorithm, based on natural selection and genetics [17]. In every GA there 
is a population of individuals, named chromosomes which a possible solution is encoded in 
each of them. A number is assigned to every chromosome showing its fitness to survive and 
reproduce. This number is calculated by fitness function. To form a new generation, chromo- 
somes are selected based on their fitness, i.e., the fitter chromosome has more chance to be 
selected for reproduction. Tournament selection is one of the selection methods which selects 
the fittest chromosome in a number of randomly selected chromosomes. Selected chromosomes 
(parents) are subjected to genetic operations of crossover and mutation to reproduce offsprings. 
Crossover is combination of two parents with a probability, Pc, creating one or two new chro- 
mosomes. Parents interchange some parts of their chromosomes at some randomly chosen 
places, e.g., in two-point crossover, two places are chosen randomly then everything between 
these places is swapped between parents. Then offsprings are passed to mutation stage which 
is a genetic operation that alters some places of each chromosome at a very low probability of 
occurrence denoted by Pm- Offsprings create a new generation and this procedure continues 
until stop condition has been reached, e.g., a solution is found or the user specified maximum 
number of generations is evolved. 



Appendix II 

For the reader convenience we present here explicit realization of generators Ai,...,A8 of Lie 
algebra su{3): 

[ I o\ 

Ai = 10 , As = i 





^ -i ^ 



y y 



^ 1 ^ 



, A3 



0-10 




FEW measure and GA 



15 



A4 = 



^001^ 





1 



A, 



/O -^^ 





1 



^ 



1 
1 



A7 = 



^ ^ 



1 

1 



A« 



1 

7! 



^10 ^ 



1 
0-2 
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Figure Captions 

Figl: Caption: E{p) versus fidelity F for Werner state p^. Parameters of program: 
Pm = 0.007, Pc = 0.7, Nga = 320, Gga = 300, A^^gjvi = 400 and Nqn2 = 5. Resolution time 
for each generation is 1 seconds on a Celeron 2 GHz PC. 

Fig2: Caption: E{p) versus q for GHZ-W mixture state pg. Parameters of program: 
Pm = 0.007, Pc = 0.7, Nga = 640, Gga = 300, Nqni = 500 and A^qjv2 = 8. Resolution time 
for each generation is 3.75 seconds on a Celeron 2 GHz PC. 

Fig3: Caption: E{p) versus a for two-qutrit state p^ . Parameters of program: P^ = 0.007, 
Pc = 0.7, Nga = 810, Gga = 300, Nqm = 800 and Nqn2 = 10. Resolution time for each 
generation is 45 seconds on a Celeron 2 GHz PC. 



